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Abstract 



We propose to use events with radiated photons in e + e collisions to measure 
I^H . the total cross section of e + e~ — ► hadrons as a function of the center of mass 

ON \ energy. The Monte Carlo simulation for the collider DAPHNE shows that a 

competitive accuracy can be achieved with this method. 

> 

ON ■ I. INTRODUCTION 

ON 

During the past years electroweak precision measurements have become one of the central 

issues in particle physics. The indirect determination of the mass of the top quark through 

ON ' its impact on quantum corrections prior to its observation in proton-antiproton collisions 

at the TEVATRON can be considered as one of the triumphs of the present theoretical 

framework. Similarly the indirect determination of the mass of the W boson is in perfect 

Ph! agreement with the present measurements at LEP and the TEVATRON, which motivates a 

c-j | combined fit to the parameters of the Standard Model. 

Recently the measurements have become sufficiently precise to even give first indications 
for a relatively small mass of the Higgs boson with an upper limit of around 200 GeV. 
Significantly improved measurements, in particularly of M\y and the left-right asymmetry 
with polarized beams at the SLC are expected to come in the next couple of years which 
might lead to a fairly precise indirect determination of the Higgs mass M#, thus repeating 
the success of the top quark mass determination. 

An important ingredient in all these fits is the fine structure constant at the mass of the 
Z boson. Its running from the Thompson limit to the high scale Mz is largely determined 
by the hadronic contribution to the vacuum polarization which in turn can be expressed via 
dispersion relation by o"had 5 the cross section for e + e~ — ► hadrons . At present the dominant 
uncertainty in the SM fits is due to the fairly large experimental error in <7h a d if the analysis is 
based on data only. Also the theoretical interpretation of future, improved determinations of 
the muon anomalous magnetic moment is affected by the limited knowledge of the hadronic 
contributions. Significant progress can only be achieved through an improved determination 
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of the cross section, either through substitution of inadequate data by theoretical predictions 
based on perturbative QCD, wherever applicable, or through improved measurements of o"h a d 
over a wide energy range, or by a combination of both. 

All present, data based evaluations make use of a large variety of experiments at different 
accelerators, a consequence of the large energy range to be spanned by the dispersion integral. 
Some of the data points, e.g. around the p resonance, are now extremely accurate, others, 
e.g. in the region between 2 and 3 GeV have large errors and are only marginally consistent 
with the predictions based on pQCD. While a highly precise scan through the whole region 
would certainly be desirable, one may, as a viable alternative, exploit the high luminosity 
of oncoming e + e~ machines which operate at fixed energy and use the radiative return to 
lower energies. This would lead to an improved determination of Chad over a large energy 
range. The feasibility of this approach is the subject of this work. 

Similar considerations can be found in |l|-||j|. However, in most of these papers only 
initial state radiation is considered which is a priori legitimate for very small angles. In fact, 
Ref. P| uses explicitly the small angle approximation and includes higher order corrections 
within the same approximation. This is an attractive option if photon detection at very 
small angles is experimentally feasible. Unfortunately, this seems to be not the case for 
the larger part of the oncoming high luminosity e + e~ machines. In the present paper we 
consider photons at somewhat larger angles, 6 1 > 7° and include initial and final state 
radiation and incorporate in addition collinear initial state radiation through the structure 
function technique. 

For the present purpose it is important to separate contributions from initial and final 
state radiation. Three techniques can be envisaged: 

i) Initial state radiation strongly dominates in events with photons emitted at small 
angles relative to the beam axis. Imposing at the same time an angular separation between 
photon and pions decreases the contribution of the final state radiation even further; 

ii) the forward-backward asymmetry of the pions which originates from the interference 
between initial and final state radiation can be measured and used to test the underlying 
model for final state radiation; 

iii) the angular distribution of photons from final state radiation involves terms pro- 
portional to sin 2 6 1 and cos 2 9 1 only and is distinctly different from the one of initial state 
radiation with its strong peaking at small angles. By defining suitable projectors it may be 
possible to distinguish the two components. 

As a particular example we consider here the case of collider DAPHNE at Frascati 
which would operate on resonance; the total energy of the collision therefore being equal 
-\/i = 1.02 GeV. The physical program of KLOE collaboration is focused on the precision 
study of the CP violation in kaon system. However, the events with radiative return to lower 
energies will be in quantity there. These events can be either considered as a background 
to an interesting physics or utilized to study the cross section e + e~ — > hadrons at lower 
energies. In what follows we show that the second option appears to be rather realistic. 
The studies presented here are based on the Monte Carlo event generator Hj that has been 
written for this purpose and which, in addition, should be useful for the analysis of the data. 









Event rates 


Collider 
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Annual luminosity, fb l 


n ro n 7° R 1 n° 

"min "J "min ' "min J-U 


DAPHNE 


1.02 


1.35 


18 • 10 6 16 ■ 10 6 14 • 10 6 


B — factory 


10.6 


100 


4-10 6 3.5 ■ 10 6 3-10 6 


LEP1 


92 


0.24 


125 109 93 


LEP2 


183 


0.2 


27 24 20 



TABLE I. Estimated number of radiative events e + e~ — » tt + tt~j for different center of mass 
energies from Eq. ([!]). The cut on the photon energy is 0.1 GeV. To calculate the annual luminosity 
we used 1 year = 10 7 sec. For LEP2 we used accumulated luminosity for the 1997 run with at 
Js = 183 GeV. 



II. BASICS OF THEORETICAL CONSIDERATION 

Before we start with a more detailed discussion of the two pion case at DAPHNE, let 
us present the order of magnitude of the event rates for a few characteristic configurations. 
The differential cross section for radiative events with # min < 6 < 180° — 9 min is given byQ: 



Cf 



da 



1 + cos 9 n 




Q 2 



cos 9 n 



(1) 



Here Q 2 is the invariant mass of the hadronic system and R(Q 2 ) is cr(e + e~ — > hadrons)/cr po i nt . 
The value of Q 2 = s — 2y/su^ is fixed by w 7 , the energy of the photon. It is important to note 
that the function in the curly brackets in Eq.(|l|) is flat for Q 2 <C s. This is similar to the 
weight function in the dispersive integral for a[Mz) and should facilitate the experimental 
determination of this important quantity. We will return to this issue below. 

As an illustrative example we consider the 7r + 7r - 7 final state, a photon angular cut 
of 5°, 7°, 10°, a realistic model for the pion form factor [0] and the typical parameters of 
four colliders: DAPHNE, 5-factory, LEP and LEP2. The respective beam energies, the 
luminosities and the event rates as obtained using Eq.fll]), are listed in Table 1. 

Another interesting class of events are those with the invariant mass of the hadronic 
system between 1.5 and 3.5 GeV. The i?-ratio in this region is very poorly known and even 
a measurement with an accuracy of 5% would be considered a significant improvement. On 
one hand, no collider will operate in this region in the near future and, on the other hand, 
the R ratio could be well determined with this method at a 5-factory. In Table 2 we list 
typical event rates, adopting for simplicity a constant R value of 2. 

It is clear from these tables that large event rates are at hand, even for severe cuts, which 
confine the photon to the region of relatively large angles. The main issues are the control 
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TABLE II. Estimated number of radiative events e + e — > hadrons + 7 at i?-factories for one 
year of running (y/s = 10.6 GeV, integrated luminosity 100 fb _ ). We use R = 2 in the interval 
1.5 < \fW < 3-5. 



of final state radiation and radiative corrections. These questions will now be studied in 
more detail for the 7r + 7r - 7 state and for the KLOE experiment. 

For yfs < 1 GeV the dominant hadronic final state produced in e + e~ annihilation consists 
of a pair of charged pions. For this reason we limit the discussion to the reaction e + e~ — > 
7r + 7r~7 and demonstrate below that by measuring the momenta of the final state particles 
at fixed center of mass energy, one can extract accurate information about the cross section 
of e + e~ — > 7r + 7r~ over a wide range of energy. 

In the reaction e + e~ — > 7r + 7r~7 the photon can be emitted either by the electron or 
positron (ISR) or by the pions in the final state (FSR). The corresponding amplitude reads: 



M = M 



ISR 



M 



FSR- 



(2) 



The calculation of ISR is a straightforward application of the Feynman rules for QED. 
In contrast, with the pions not being elementary, the calculation of FSR is in general more 
tricky. In the current version of our Monte Carlo program, we consider pions as point-like 
particles as far as FSR is concerned and use the standard Feynman rules of scalar QED. 
Being definitely not rigorous, this approach should provide a reasonable estimate of FSR for 
photons emitted at small angles relative to pions and an upper limit for photon emission 
at large angles. Moreover, below we present a useful tool to control the accuracy of this 
approximation. Final state radiation through the radiative decay <fi — > 7r + 7r~7 is not included 
in the analysis (for a recent discussion see and references therein). It can be considered 
as a special case of FSR and we expect that this effect can also be controlled through its 
interference with the dominant amplitude. 

The interaction of the virtual photon with the 7r + 7r~ system is described by a pion form 
factor F-xiQ 2 ), where Q 2 is the virtuality of the photon. A convenient parameterization of 
this form factor was proposed in Ref . . The pion form factor is described by a sum of the 
Breit-Wigner resonances with a Q 2 dependent width: 



FAQ 2 ) 



BWp (l+«BW„) + j3BWp/ + lBWf) 

1 + /3 + 7 



2 
Tfl 

BW P (Q 2 ) = 2 n2 P f7uv jm - (3) 

m 2 p -Q 2 - iy/Q 2 T p (Q 2 ) 

In the above equation a, /3, 7 are the parameters of the model (see Ref. M). For Q 2 = 
one obtains -^(0) = 1; that ensures the proper charge of a pion. In what follows we shall 
use this parameterization. 

In general, by measuring the photon energy uj^ one determines the invariant mass of the 
pions. Consider now the square of the ISR amplitude. In this case the differential cross 
section is evidently proportional to the pion form factor squared at the proper momentum 
scale: 

(■*?») ~ l F ^ 2 )l 2 ~ <w^+.-(Q 2 ). (4) 

\ a( ^ /ISR 

A different situation occurs when the photon is emitted from the final state. In this case 
the differential cross section is proportional to 

/ dn * ~|^)lV0-e+e-^+.-(Q 2 ). (5) 



dQ 



FSR 



Clearly, in this case measuring the energy of the final state photon would not help in 
determining the proper energy of the collision and therefore the cross section of e + e~ — > 
7r + 7r~. For this reason FSR and ISR/FSR interference should be regarded as a "background" 
and the contribution due to ISR as a "signal". It is definitely not possible to "switch 
off" FSR completely , hence FSR and the interference will lead to a systematic error in 
the determination of a(e + e~ — > 7r + 7r~) from a(e + e~ — > ir + 7t~ , y). The choice of a proper 
selection criterion is thus the only possibility to suppress the contribution of FSR. In fact, 
the DAPHNE luminosity is so high, that even after imposing quite severe cuts, the number 
of radiative events remains huge, rendering the statistical accuracy of the measurement a 
minor problem. For this reason, it is the level of suppression of FSR and the control of the 
remainder that determines the accuracy of the measurement. For the tt + tt~ final state, recent 
measurements in Novosibirsk have pushed the uncertainty in the pion form factor down to 
one per cent for ^/s < 1 GeV |§. Hence, for the proposed method to be competitive, the 
level of the FSR contribution should not exceed this number. 

In addition the FSR amplitude can be calibrated experimentally by measuring the ISR- 
FSR interference which is odd under charge conjugation and thus gives rise to a relatively 
large forward-backward asymmetry. 

In practice, the situation is further complicated by additional collinear initial state radia- 
tion which changes the total energy of the collision and can not be detected. This effect must 
be taken into account in any realistic Monte Carlo event generator. As usual, the collinear 
initial state radiation is described by the structure functions which provide a resummation of 
the logarithmic corrections O (a log 2 (s/ m 2 e )) to all orders in the coupling constant. We also 
assume that the invariant mass of the 7r + 7r~7 system can be measured sufficiently accurate. 
Events with the invariant mass of the final state much smaller than the total energy of the 
collision squared can thus be rejected without loosing too much luminosity. This allows to 
reduce significantly the kinematic distortion of the events due to collinear initial state radi- 
ation. In the Monte Carlo program this is realized by requiring a minimal invariant mass of 
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the final state 7T + 7r _ 7. For the structure functions the formulas of Ref. || are used. For the 
purpose of illustration, we present the results for the pion invariant mass distribution with 
and without collinear initial state radiation in Fig.[l]. 

In what follows we discuss two issues essential for a successful measurement. First, we 
show that by choosing a proper cuts one can significantly suppress the contribution due to 
FSR. Second, we demonstrate that by measuring the forward-backward asymmetry of the 
produced pions one can check how realistic the final state radiation amplitude is described. 
Confronting our predictions for the forward-backward asymmetry with the measurements, 
one can gain confidence in the predictions based on our Monte Carlo event generator. If 
not stated otherwise, the additional collinear radiation is incorporated with the structure 
function technique. 

III. ISR DOMINANCE 

The spectrum of photons emitted by an ultrarelativistic particle is described by the 
well-known formula: 

dcr oc d<7 - 



UJ.., 



v, 



Hence the relativistic particle strongly radiates in a small angular cone along its direction of 
motion. For this reason, a suppression of the FSR is achieved by selecting events where the 
pions are geometrically separated from the hard photon. Simultaneously, an enhancement 
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of the ISR is achieved if one selects the events where the photon is emitted at small angles 
relative to the collision axis. In reality, the KLOE detector seems to be able to measure 
photons down to 9 1 ~ 7°. Therefore, an appropriate configuration is, for example, obtained 
by requiring 30° < 9 W < 150° for pions and 7° < # 7 < 20° or 160° < 6 1 < 173° for photons. 
The prediction for the Q\ +1X ^ invariant mass distribution is shown in Fig.|2|. Note that 
we have included events with the photons emitted both in the forward and the backward 
direction. Nevertheless, a forward backward asymmetry of the tt + (or, with opposite sign, of 
the 7r~) will be observed, which should not be confused with the trivial kinematic asymmetry 
if we would consider photons with 7° < # 7 < 20° only. 

It follows from Fig.0 that for the chosen cuts the contribution of the FSR radiation to 
the total cross section is of the order of 1 per cent only. Also, one does not loose too much 
in the event rate since the pion angular distribution is relatively uniform and the photon 
angular distribution is strongly peaked at small angles. 

Note that the contribution of the FSR to such a configuration comes only from the 
emission of photons by pions at large angles. In this case the pions in the intermediate 
state are far off-shell. For this reason, one expects that the point-like interaction of pions 
to photons is not likely to be an adequate description. However, we believe, that in this 
case the point-like interaction overestimates the strength of the interaction and therefore 
gives a larger contribution due to the FSR. Still, it is desirable to check how reasonable the 
interaction of photons to pions is described by our ansatz. A possible check is considered in 
the next Section. 

The distribution shown in Fig.|^ may be compared with the distribution when no geo- 
metric separation is imposed on photons and pions FigJ^. In this case, one finds a significant 
FSR contribution. An interesting fact to note, is that the FSR contribution decreases for 
smaller Ql +7T - ■ This implies that possibly better results can be achieved if a combined Q\+^- 
- cos# 7 cut is applied. Large Q 2 , on the other hand, correspond to soft photons where our 
description of FSR is expected to work reasonably well. 

Let us imagine that the application of the proper selection criteria resulted in a suppres- 
sion of the final state radiation. In this case, even measuring the total cross section integrated 
over the invariant masses of pions could provide useful information for the determination of 
the fine structure constant. 

To illustrate this point, we consider the expression for the cross section due to ISR: 

<*—> _ |V (Q) |s u _ jgf i (*; + <w i , - i^gn . ( 6) 



dQ 2 dcosfl 7 3s 2 ' nK ^ n \ Q 2 J |Q 2 (s-Q 2 )l-cos 2 # 7 

Consider now the case s ~^> Q 2 , i.e. a creation of the pion pair with the invariant mass 
much smaller than the energy of the collision. The above equation is then simplified: 

da ee ^ W7T1 a (l + cos# 2 ) 



dQ 2 dcosfl 7 2irs "^ '(l-cos6 2 ) 

The total cross section for e + e~ — > 7r + 7r _ 7 is obtained if one integrates over the invariant 
mass of pions and the photon production angle: 




FIG. 2. A comparison of the ISR contribution (dashed line) with the complete result (solid 
line). Plotted is the cross section da(e + e~ — ► 7r + 7r~7)/dQ 2 + _ in nbarn/GeV 2 as a function 
of Ql+ W - in GeV 2 for ^ = 1.02 GeV. The cuts are 7° < 7 < 20° or 160° < 6> 7 < 183°, 
30° < 9^ < 150°, w 7 > 0.02 GeV and the invariant mass of the detected particles in the final state 



Q 2 + _ > 0.9 GeV 2 . 
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FIG. 3. The same as in Fig.|, but for angular cuts 30° < 6~, 9* < 150°. 
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FIG. 4. Pion (vr + ) angular distribution in the laboratory frame. The cut on the photon angle 
is 60° < 6 y < 120°. 
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This expression should be compared with the small s contribution to fthad- Considering 
the contribution of the low-energy region, s <C M§, one obtains: 
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We therefore conclude, that even the integrated cross section for e + e~ — > hadrons + 7 
can be effectively used to obtain the contribution of the low energy region to the value of 
the fine structure constant at the Z resonance, provided that the suppression of the final 
state radiation is achieved using appropriate selection criteria. 

In a realistic situation, the accurate determination of <r(e + e~ — > 7r + 7r~) from a(e + e~ — > 
7r + 7r~7) requires the knowledge of the contribution due to FSR and the effect of the collinear 
ISR. As we have explained above, FSR can be significantly suppressed using appropriate 
cuts on the photon and pion emission angles. After this is achieved, an appropriate strategy 
which can be used to extract the pion form factor would be to fit the parameters of the 
model for the pion form factor described in [[?] to an observed Q\+^- -distribution using our 
Monte Carlo event generator which incorporates all the effects related to the collinear initial 
state radiation. 
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FIG. 5. Pion (vr + ) angular distribution in the laboratory frame. The photon angle is restricted 
to 7° < 7 < 20° or 160° < 7 < 183°. The cut on the pion angle is 30° < 0„.+ )7r - < 150°. 



IV. TESTING THE MODEL FOR THE FSR 

The cross section for e + e~ — > 7r + 7r~7 can be written as: 

da k\M\ 2 = |Msr| 2 + I-MfsrJ 2 + 2Re [MsrA*£ S r] • (9) 

We have explicitly separated the contribution due to ISR, FSR and the interference of the 
two amplitudes. 

If the photon is emitted from the initial(fmal) state, the pion pair is produced with charge 
parity C = — 1(+1). Therefore, the contribution of the third term in Eq.(|9|) vanishes if one 
integrates over the kinematic variables of the pions. However, this term causes a significant 
charge asymmetry and, correspondently, a forward-backward asymmetry of tc + and 7r~. 

As usual, the asymmetry is defined as: 
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J a cos _i_ n J a cos & 4. " 



(10) 



The asymmetry is also a function of the cuts on the photon angle and energy which is 
implicitly assumed in the above equation. This asymmetry should not be confused with 
the trivial asymmetry from kinematics which arises if one selects photons in one hemisphere 
only. 

For symmetric cuts on the photon angular distribution, the numerator of this expression 
is non zero only due to the interference term in Eq.(^J), the denominator equals to the total 
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cross section. Given the dominance of the initial state radiation in the total cross section, the 
forward-backward asymmetry is (roughly) proportional to the ratio .MfsrV-Misr an d thus 
is a direct measure of the final state radiation. For this reason, measuring the asymmetry 
would give an effective control over the accuracy of the ansatz used to describe the interaction 
of photons to pions. 

To demonstrate the magnitude of the effect consider the situation, where the standard 
cuts on the photon energy and the invariant mass of 7r + 7r~7 are applied. The photon 
angular cuts are such that 60° < 9 1 < 120°, also lo 1 > 0.02 GeV. The angular distribution 
for positively charged pions is then presented in Fig.|j. One sees, that the distribution is not 
symmetric: positively charged pions are mainly emitted in the backward direction. Using 
Eq . (|10"1) , we then obtain that the asymmetry A equals 20%. The total cross section in this 
case equals 3.29 nbarn which implies that up to 3.29 • 10 6 events will be observed in a year 
of running. Without any doubt, the twenty percent effect will be clearly visible with such a 
statistic. 

Alternatively, we present in Fig.|5] the angular distribution of pions for the case where 
the photon is emitted at small angles, well separated from pions: 7° < 6* 7 < 20° or 160° < 
6> 7 < 183°, 30° < 9 n + tn - < 150° and co> 7 > 0.02 GeV. In this case the asymmetry is reduced 
to A = 1.5%, demonstrating again the smallness of final state radiation for this case. 

The separation between initial and final state radiation is also possible for the symmetric 
piece on the basis of the distinctly different angular distributions. Let us for the moment 
ignore the folding with the radiator functions and consider strictly the reaction e + e~ — ► 
7r + 7r~7. Final state radiation alone leads necessarily to a photon angular distribution of 
the form a + /3cos 2 # 7 , initial state radiation on the other hand is described by Eq.(^). 
Given sufficient statistics, it is clear that a combined fit will allow the separation of the two 
components. 

V. CONCLUSIONS 

In this paper we have argued, that one can effectively use the radiative return to low 
energies in order to measure the cross section of e + e~ — > hadrons at variable center of mass 
energy, of course lower than the energy of the colliding beams. We have studied the feasibility 
of this approach with particular emphasis on the two pion final state - a typical situation 
for the DAPHNE ring. We have produced a realistic tool for these studies, a Monte Carlo 
event generator, which provides an option for choosing the cuts and includes both initial, 
final and collinear initial radiation. 

By using special selection criteria and measuring the energies and momenta of the par- 
ticles in the final state one can suppress the contribution of the final state radiation to the 
one per cent level. We have also demonstrated that there will be a sizable forward-backward 
asymmetry, which can be used to check the quality of our description of the final state 
radiation. 

For the forward-backward asymmetry to be significant and to allow for this check, one 
should keep the photon angle sufficiently large, otherwise the ratio A^isr/A^fsr becomes too 
small. The opposite requirement must be fulfilled for the accurate extraction of a(e + e~ — > 
hadrons), as discussed in Sect. 3. For this reason, all events with photons appear to be useful 

11 



for extracting interesting information. 

The proposed technique was mainly discussed in connection with the DAPHNE collider. 
We would like to stress, however, that the use of the radiation return to lower energies can 
be used in a more general context and the same technique should, in principle, be applicable 
also for other "particle factories" . 
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